Microscopic analysis of sub-barrier fusion enhancement in 132 Sn+ 40 Ca vs. 132 Sn+ 48 Ca 
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We provide a theoretical analysis of recently measured fusion cross sections which show a sur- 
prising enhancement at low E c , m , energies for the system 132 Sn+ 40 Ca as compared to the more 
neutron-rich system 132 Sn+ 48 Ca. Dynamic microscopic calculations are carried out on a three- 
dimensional lattice with a time-dependent density-constrained density functional theory. There are 
no adjustable parameters, the only input is the Skyrme effective NN interaction. Heavy-ion poten- 
tials V(R), coordinate-dependent mass parameters M(R), and total fusion cross sections a(E c . Ta .) 
are calculated for both systems. We are able to explain the measured fusion enhancement in terms 
of the narrower width of the ion-ion potential for 132 Sn+ 40 Ca, while the barrier heights and positions 
are approximately the same in both systems. 
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I. INTRODUCTION 

Radioactive ion beam facilities enable us to study fu- 
sion reactions with exotic neutron-rich nuclei. An im- 
portant goal of these experiments is to study the effects 
of neutron excess (N — Z) on fusion. In several exper- 
iments, large enhancements of sub-barrier fusion yields 
have been observed for systems with positive Q values 
for neutron transfer. Recently, at the HRIBF facility a 
series of experiments has been carried out with radioac- 
tive 132 Sn beams and with stable 124 Sn beams on 40 ' 48 Ca 
targets PQ. It turns out that the 40 Ca+Sn systems have 
many positive Q values for neutron-pickup while all the 
Q values for 48 Ca+Sn are negative. However, the data 
analysis reveals that the fusion enhancement is not pro- 
portional to the magnitudes of those Q values. 

Particularly puzzling is the experimental observa- 
tion of a sub-barrier fusion enhancement in the system 
132 Sn+ 40 Ca as compared to more neutron-rich system 
132 Sn+ 48 Ca. This is difficult to understand because the 
8 additional neutrons in 48 Ca should increase the attrac- 
tive strong nuclear interaction and thus lower the fusion 
barrier, resulting in an enhanced sub-barrier fusion cross 
section. But the opposite is found experimentally. A 
coupled channel analysis [T] of the fusion data with phc- 
nomenological heavy-ion potentials yields cross sections 
that are one order of magnitude too small at sub-barrier 
energies, despite the fact that these ion-ion potentials 
contain many adjustable parameters. 

The time-dependent Hartree-Fock (TDHF) theory pro- 
vides a useful foundation for a fully microscopic many- 
body theory of large amplitude collective motion [21 13| 
including deep-inelastic and fusion reactions. But only 
in recent years has it become feasible to perform TDHF 
calculations on a 3D Cartesian grid without any symme- 
try restrictions and with much more accurate numerical 
methods [3HH]- In addition, the quality of effective inter- 
actions has been substantially improved [5lU2|. During 
the past several years, we have developed the Density 
Constrained Time-Dependent Hartree-Fock (DC-TDHF) 
method for calculating heavy-ion potentials |13j . and we 



have applied this method to calculate fusion and capture 
cross sections above and below the barrier. As of to date, 
we have studied a total of 18 different systems, including 

132,124 Sn+ 96 Zr 132 Sn+ 64 Ni flfiftft 16 O+ 208p b ^ 

and hot and cold fusion reactions leading to superheavy 
element Z = 112 [TS]. Most recently, we have investi- 
gated sub-barrier fusion and pre-equilibrium giant res- 
onance excitation between various calcium + calcium 
isotopes [15], and between isotopes of oxygen and car- 
bon [2DJ that occur in the neutron star crust. In all cases, 
we have found good agreement between the measured fu- 
sion cross sections and the DC-TDHF results. This is 
rather remarkable given the fact that the only input in 
DC-TDHF is the Skyrme effective N-N interaction, and 
there are no adjustable parameters. 

This paper is organized as follows: in Section II we 
summarize the theoretical formalism and show results 
for the ion-ion potentials calculated with the DC-TDHF 
method. In Section III, we discuss the corresponding 
total fusion cross sections, and we explain the observed 
fusion enhancement in terms of the narrower width of 
the ion-ion potential for 132 Sn-(- 40 Ca. Our conclusions 
are presented in Section IV. 



II. THEORETICAL FORMALISM AND 
ION-ION POTENTIALS 

Currently, a true quantum many-body theory of bar- 
rier tunneling does not exist. All sub-barrier fusion calcu- 
lations assume that there exists an ion-ion potential V(R) 
which depends on the internuclear distance R. Most of 
the theoretical fusion studies are carried out with the 
coupled-channels (CC) method PTH2"1] in which one uses 
empirical ion-ion potentials (typically Woods-Saxon po- 
tentials, or double-folding potentials with frozen nuclear 
densities). 

While phenomenological methods provide a useful 
starting point for the analysis of fusion data, it is de- 
sirable to use a quantum many-body approach which 
properly describes the underlying nuclear shell structure 
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of the reaction system. We have developed a micro- 
scopic approach to extract heavy-ion interaction poten- 
tials V(R) from the TDHF time-evolution of the dinu- 
clear system. The interaction potentials calculated with 
the DC-TDHF method incorporate all of the dynamical 
entrance channel effects such as neck formation, parti- 
cle exchange, internal excitations, and deformation ef- 
fects [35] . While the outer part of the potential barrier 
is largely determined by the entrance channel properties, 
the inner part of the potential barrier is strongly sen- 
sitive to dynamical effects such as particle transfer and 
neck formation. 

The TDHF equations for the single-particle wave func- 
tions 

fc({^.})^A(r,f) = ift~fo(r,i) {X = 1,...,A) , (1) 

can be derived from a variational principle. The main ap- 
proximation in TDHF is that the many-body wave func- 
tion $(i) is assumed to be a single time-dependent Slater 
determinant which consists of an anti-symmetrized prod- 
uct of single-particle wave functions 

$(n,...,r A ;t) = (Al)- 1 / 2 det\<j> x {n,t)\ . (2) 

In the present TDHF calculations we use the Skyrme 
SLy4 interaction [9] for the nucleons including all of the 
time-odd terms in the mean-held Hamiltonian [5 . The 
numerical calculations are carried out on a 3D Cartesian 
lattice. For 40 ' 48 Ca+ 132 Sn the lattice spans 50 fm along 
the collision axis and 30 — 42 fm in the other two di- 
rections, depending on the impact parameter. First we 
generate very accurate static HF wave functions for the 
two nuclei on the 3D grid. In the second step, we apply a 
boost operator to the single-particle wave functions. The 
time-propagation is carried out using a Taylor series ex- 
pansion (up to orders 10 — 12) of the unitary mean-field 
propagator, with a time step At = 0.4 fm/c. 

In our DC-TDHF approach, the time-evolution takes 
place with no restrictions. At certain times t or, equiva- 
lently, at certain internuclear distances R(t) the instan- 
taneous TDHF density 

PtdhfM) =< $(t)|p|$(*) > (3) 

is used to perform a static Hartree-Fock energy minimiza- 
tion 

5 <<& P \H- J d 3 r A(r) p(r) |$ p >= (4) 

while constraining the proton and neutron densities to 
be equal to the instantaneous TDHF densities 

< $p\p\$ p >= Ptdhf(j, t) . (5) 

These equations determine the state vector $ p . This 
means we allow the single-particle wave functions to re- 
arrange themselves in such a way that the total energy 
is minimized, subject to the TDHF density constraint. 



In a typical DC-TDHF run, we utilize a few thousand 
time steps, and the density constraint is applied every 
10 — 20 time steps. We refer to the minimized energy as 
the "density constrained energy" E^c{R) 

Evc(R) =< ^ P \H\^ P > . (6) 

The ion-ion interaction potential V(R) is essentially the 
same as Edc{R), except that it is renormalized by sub- 
tracting the constant binding energies E Al and Ej± 2 of 
the two individual nuclei 

V(R) = E DC {R) -E Al -E A2 . (7) 

In Fig. [I] we compare the heavy-ion interaction poten- 
tials V(R) for the systems 132 Sn+ 40 > 48 Ca. It should be 
noted that DC-TDHF contains particle transfer "on av- 
erage" [12], but it does not describe individual transfer 
channels. We find the unexpected result that the barrier 
heights and positions are approximately the same in both 
cases, but the barrier width for 132 Sn+ 40 Ca is substan- 
tially smaller. 




R(fm) 

FIG. 1. DC-TDHF calculation of the heavy-ion potentials 
V(R) for the systems 132 Sn+ 40 ' 48 Ca. The ion-ion potentials 
are energy-dependent and have been calculated at -Etdhf = 
120 MeV. 

Using TDHF dynamics, it is also possible to compute 
the corresponding coordinate dependent mass parame- 
ter M(R) [17] , At large distance R, the mass M(R) is 
equal to the reduced mass y, of the system. At smaller 
distances, when the nuclei overlap, the mass parame- 
ter generally increases. We find that the structure of 
M (R) for the 132 Sn+ 40 Ca reaction is fairly similar to the 
mass parameter calculated for 132 Sn+ 48 Ca (see Fig. 6 of 
Ref. [35]), and it is therefore not shown here. 

Instead of solving the Schrodinger equation with coor- 
dinate dependent mass parameter M(R) for the heavy- 
ion potential V(R), it is numerically advantageous to 
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use the constant reduced mass [i and to transfer the 
coordinate-dependence of the mass to a scaled potential 
U (R) using a scale transformation 



(R, M(R), V(R)) — > (R, [i, U(R)) . 



(8) 



Details are given in Ref. [T7] . 

In Fig. |2] we display the transformed potentials U(R) 
which correspond to the constant reduced mass /1. 
A comparison of Fig. [T] and Fig. [2] reveals that the 
coordinate-dependent mass changes only the interior re- 
gion of the potential barriers, and this change is most 
pronounced at low E c . m . energies. Note that the trans- 
formation to a constant mass parameter preserves the 
basic features found for the original potentials, i.e. the 
narrower width of the ion-ion potential for Sn+ Ca, 
while the barrier heights and positions are approximately 
the same in both systems. 
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FIG. 2. Transformed heavy-ion potentials U(R) correspond- 
ing to the constant reduced mass {j,. 



at low energy (-Etdhf = 120 MeV) has a barrier of only 
E B = 112.3 MeV located at R = 12.4 fm. 
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FIG. 3. (Color online) Solid lines: original heavy-ion poten- 
tials V(R). Dashed lines: transformed potentials U(R) cor- 
responding to the reduced mass /j,. The potentials have been 
calculated at three different energies. 

The Schrodinger equation corresponding to the con- 
stant reduced mass /i and the scaled potential U (R) has 
the familiar form 



-h 2 d 2 



h 2 e(e + 1) 



U{R) - E c 



MR) = 



2[i dR 2 2/ii? 2 

(9) 

By numerical integration of Eq. Q using the 
well-established Incoming Wave Boundary Condition 
(IWBC) method [ST] we obtain the barrier penetrabil- 
ities To which determine the total fusion cross section 



CfusC^c 
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2/J,E c 



J2(^ + l)Te(E c 



(10) 
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In general, our DC-TDHF calculations show that ion- 
ion potentials for heavy systems are strongly energy- 
dependent. By contrast, DC-TDHF calculations for light 
ion systems such as 16 0+ 16 show almost no energy- 
dependence even if we increase E c . m _ by a factor of 
four |27) . Even in reactions between a light and a very 
heavy nucleus such as 16 O+ 208 Pb, we see only a rela- 
tively weak energy dependence of the barrier height and 
width 17 . In Fig. [3] the original potentials V(R) (solid 
lines) and the transformed potentials U(R) are shown at 
three different TDHF energies. We notice that in these 
heavy systems the potential barrier height increases dra- 
matically with increasing energy Etbhf, and the bar- 
rier peak moves inward towards smaller internuclear dis- 
tances. The potential U(R) calculated at high energy 
(Stdhf =_ 180 MeV) has a barrier E B = 115.3 MeV lo- 
cated at R = 11.8 fm, whereas the potential calculated 



III. FUSION CROSS SECTIONS 

In Figures [4] and [5] we show fusion cross sections mea- 
sured at HRIBF [JJ for the systems 132 Sn+ 40 < 48 Ca. A 
comparison of the fusion cross sections at low energies 
yields the surprising result that fusion of 132 Sn with 40 Ca 
yields a larger cross section than with 48 Ca. For exam- 
ple, at E c , m , = 110 MeV we find an experimental cross 
section of w 6 mb for 132 Sn+ 40 Ca as compared to 0.8 
mb for the more neutron-rich system 132 Sn+ Ca. If the 
data are scaled for trivial size effects (nuclear radii) the 
difference between the "reduced" cross sections is found 
to be even larger, see Fig. 6 of Ref. [T]. The experimen- 
talists have carried out a coupled channel (CC) analysis 
of the fusion data with phenomenological Woods-Saxon 
potentials which generally underpredict the data at low 
E c .m. energies. In the case of 132 Sn+ 40 Ca, the CC model 
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calculations yield cross sections which differ by a fac- 
tor of 10 or more from the data, despite the fact that 
the optical model potentials contain 7 adjustable pa- 
rameters. Interestingly, it is possible to get a good fit 
to the data using the empirical Wong model (tunneling 
through a single parabolic barrier, with 3 adjustable pa- 
rameters). In this case the analysis reveals an unusually 
large curvature of the barrier, Hui = 13.13±1.09 MeV, for 
132 Sn+ 40 Ca as compared to the 132 Sn+ 48 Ca system with 
only Huj = 5.77 ± 0.52 MeV. A large curvature implies a 
narrow parabolic barrier. Of course, these values are sim- 
ply fits to the measured data; the model does not explain 
why the barrier curvatures are so dramatically different. 
As we will see, our microscopic DC-TDHF theory de- 
scribes these results naturally in terms of the underlying 
mean-field dynamics, without any adjustable parameters. 
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FIG. 4. (Color online) Total fusion cross sections for 
132 Sn4- 40 Ca. The cross sections have been calculated with 
the DC-TDHF method for several energy-dependent ion-ion 
potentials U (R) some of which are displayed in Fig. [3] The 
experimental data are taken from Ref. [T]. 

In Fig. [4] we show the excitation function of the total 
fusion cross section for 132 Sn+ 40 Ca. The cross sections 
have been calculated with five different energy-dependent 
ion-ion potentials; the corresponding DC-TDHF energies 
are listed in the figure. The main point of the chosen 
display is to demonstrate that the energy-dependence of 
the heavy-ion potential is crucial for an understanding 
of the strong fusion enhancement at subbarrier energies: 
At very high energy (-Etdhf = 180 MeV) the potential 
approaches the limit of the frozen density approximation: 
the collision is so fast that the nuclei have no time to 
rearrange their densities. We observe that the measured 
fusion cross sections at energies E c , m , > 118 MeV are 
well-described by this high-energy potential. 

As we have seen in Fig. [3j the potential barrier height 
decreases dramatically as we lower the energy -Etdhf- 
Due to the slow motion of the nuclei at sub-barrier en- 
ergies, the nuclear densities have time to rearrange re- 
sulting in neck formation, neutron transfer, and surface 



vibrations. All of these effects are included in DC-TDHF, 
and apparently they reduce the fusion barrier height and 
strongly modify the interior region of the heavy-ion inter- 
action potential. As a result of the decrease in the fusion 
barrier, one finds strongly enhanced fusion cross sections 
at low sub-barrier energies. In particular, we observe 
that the data points measured at the two lowest energies 
_E c m = 108.6 and 111 MeV are described well by the 
heavy- ion potential calculated at Etdhf = 115 MeV. 
This is the lowest-energy potential U(R) we have been 
able to calculate using the DC-TDHF method, with a 
potential barrier Eb = 111.5 MeV located at R = 12.4 
fm. 
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FIG. 5. (Color online) Total fusion cross section obtained 
with the DC-TDHF method for 132 Sn+ 48 Ca. Cross sections 
calculated for several energy-dependent ion-ion potentials [26| 
have been interpolated in this case (single blue line). The 
experimental data are taken from Ref. [Tj. 



In Fig. [5] we show total fusion cross sections for 
132 Sn+ 48 Ca which contains 8 additional neutrons. In this 
case, we have interpolated the theoretical cross sections 
obtained with the energy-dependent DC-TDHF poten- 
tials We can see that our theoretical cross sections 
agree remarkably well with the experimental data. The 
main experimental puzzle, i.e. the fact that the low- 
energy sub-barrier fusion cross section for 132 Sn+ 40 Ca is 
substantially enhanced as compared to the more neutron- 
rich system 132 Sn+ 48 Ca, can be understood by examin- 
ing the transformed ion-ion potential shown in Fig. [2] 
Both systems are found to have approximately the same 
barrier heights and positions, but the barrier for 40 Ca has 
a narrower width, resulting in enhanced fusion. Our DC- 
TDHF approach naturally explains the results of the ex- 
perimental data analysis which used the empirical Wong 
model fit (see remarks at the beginning of this section). 
However, the microscopic potential barrier is not a simple 
parabola, and the DC-TDHF ion-ion potential is found to 
depend strongly on the energy £ c , m , for heavy systems. 
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IV. SUMMARY 



In this paper, we calculate heavy-ion interaction po- 
tentials and total fusion cross sections for 132 Sn+ 40 Ca 
at energies E c m below and above the Coulomb barrier, 
and we compare these results to the more neutron-rich 
system 132 Sn+ 48 Ca studied earlier (35]. The ion- ion po- 
tential calculations are carried out utilizing a dynamic 
microscopic approach, the Density Constrained Time- 
Dependent Hartree-Fock (DC-TDHF) method. The 
single-particle wave functions are generated on a 3D 
Cartesian lattice which spans 50 fm along the collision 
axis and 30 — 42 fm in the other two directions. The 
only input is the Skyrme N-N interaction, there are no 
adjustable parameters. 

The main objective of this paper is to give a theoreti- 
cal analysis of fusion cross sections which were measured 
recently at HRIBF. The experimental data show a sur- 
prising enhancement at low -E c .m. energies for the system 
132 Sn+ 40 Ca as compared to the more neutron-rich sys- 
tem 132 Sn+ 48 Ca. Based on geometric considerations, one 
would expect the opposite: as a result of the increased 
nuclear radius for 48 Ca, the fusion barrier for 132 Sn+ 48 Ca 
should be reduced which in turn should increase the fu- 
sion cross section. Using the microscopic DC-TDHF ap- 
proach we are able to explain the measured sub-barrier 



fusion enhancement in terms of the narrower width of 
the ion-ion potential for 132 Sn+ 40 Ca, while the barrier 
heights and positions are approximately the same in both 
systems. 

While for the fusion of light nuclei the microscopic ion- 
ion potentials are almost independent of the cm. energy, 
for heavier systems a strong energy-dependence is ob- 
served. With increasing cm. energy, the height of the 
potential barrier increases, and the barrier peak moves in- 
ward towards smaller internuclear distances (see Fig. |3| . 
This behavior of the ion-ion potential has a dramatic 
influence on the sub-barrier fusion cross sections. For 
the system 132 Sn+ 40 Ca, we have calculated heavy-ion in- 
teraction potentials at 8 different energies ranging from 
E TDHF = 115 MeV to S TD HF = 180 MeV. The time- 
dependent and density constraint calculations are com- 
putationally very intensive. The total CPU time in this 
case amounts to 192 days on a single Intel Xeon proces- 
sor. Our calculations are performed on a Dell LINUX 
workstation with 12 processors using OpenMP, which re- 
duces this time to about 16 days. 
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